J/ip production in direct and resolved 77 collisions 
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The production of J/ip mesons in 77 collisions allows for a test of factorization in color-octet processes as 
predicted by NRQCD and observed at the Tevatron. We calculate the cross sections for J/ip production with 
direct and resolved photons, including also the feed-down from \ c j and ip' decays. Our NRQCD predictions are 
nicely confirmed by recent data from the DELPHI collaboration at CERN LEP2. 



1. MOTIVATION 

One of the most interesting features of quan- 
tum chromodynamics (QCD) is the apparent dis- 
crepancy between the fundamental role of color as 
a conserved quantum number carried by quarks 
and gluons and the non-observation of color in 
physical hadronic states. As a non-perturbative 
phenomenon, the confinement mechanism re- 
sponsible for the color-singlet nature of hadrons 
is still not fully understood. One can hope to gain 
more insight by studying the relatively simple 
bound states of heavy charm and bottom quarks, 
such as J/ip or T mesons. 

Ever since QCD emerged as the fundamental 
theory of the strong interaction, models have been 
designed to explain the compensation of color in 
heavy quarkonia. Well-known examples are the 
Color Evaporation jJjlB || and Hard Comover 
Scattering Models 3' which assume that color 
is restored by soft and hard gluon exchange with 
the underlying event, respectively, and the Color 
Singlet Model (CSM) [ @, 0, which requires the 
heavy quark pair to be produced in a color singlet 
state. The CSM leads to strong selection rules 
and connects the production cross section with 
the quarkonium wave function. Today, a rigorous 
effective field theory exists in the form of non- 
relativistic QCD (NRQCD) [ |j, which utilizes a 
double expansion in the strong coupling constant 
a s and the relative quark-antiquark velocity v in 
order to factorize the hard production from the 
soft binding process [ Left-over singularities 
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in the CSM can be removed systematically into 
non-perturbative color-octet operator matrix el- 
ements (OMEs), and their numerical values can 
be fitted to describe the hadroproduction cross 
sections observed at the Tevatron. However, it 
is necessary to demonstrate the universality of 
these OMEs in other production processes, such 
as 77 collisions. Since the theoretical uncertain- 
ties from scale variations in leading order (LO) of 
a s arc quite substantial, it is furthermore impor- 
tant to include also contributions from next-to- 
leading order (NLO) virtual loop and real emis- 
sion processes. In this paper, we report on recent 
progress along these lines. 

2. REAL CORRECTIONS TO J/ip PRO- 
DUCTION IN DIRECT 77 COLLI- 
SIONS 

In direct 77 collisions, J/ip mesons with mass 
M — 2m c and finite transverse momentum px 
are produced in association with either a photon 
(see Fig. P or a jet, where the final state photon 
in Fig. [l] has to be replaced by a gluon. In the 
first case, the physical color-singlet state ( 3 5'[ 1 ' 
in spectroscopic notation) is produced directly, 
whereas in the second case the intermediate color- 
octet state ( 3 >5[ 8 ') transforms into the physical 
J /ill through soft gluon emission. Although the 
strong coupling of the hard gluon enhances the 
second process, it is strongly suppressed by the 
relevant color-octet OME, which is subleading in 
v. In NLO of a s , an unresolved dijet system with 
invariant mass Sjj, originating from two gluons 
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Figure 1. Feynman diagrams pertinent to the partonic subprocess 77 — ► J/V'7- 



or a light quark- antiquark pair, allows for the 
production of intermediate 1 S^\ 3 sf\ and 3 Pj 8 ' 
states. At small pr, virtual loop corrections have 
to be included to cancel the unphysical depen- 
dence on the cut-off Sjj > M 2 , but at large pr 
these contributions become unimportant. As can 
be observed clearly in Fig. ||, the cut-off depen- 
dence is reduced and the corrections from real 
particle emission become large [ |l(|. At small 
Pt, one does, however, not expect large radiative 
corrections, as has also been shown in a full NLO 
calculation for direct color-singlet Jji\> photopro- 
duction [H. 



3. COMPARISON OF J/i/} PRODUC- 
TION IN COMPLETE 77 COLLI- 
SIONS WITH CERN LEP2 DATA 

J/ip production in 77 collisions proceeds not 
only through direct interactions of photons with 
charm quarks, but also through resolved pro- 
cesses, where one or two of the photons emit 
quarks and gluons which then interact with the 
heavy quarks. In the small-p^ range recently 
studied by the DELPHI experiment at CERN 
LEP2 [ [12]], J/tp production is dominated by the 
single-resolved process 75 — > J/ipg, where one ini- 
tial state photon and the final state photon in Fig. 
[l] have to be replaced with gluons. In addition, 
processes with non-abelian coupling of the glu- 
ons and processes involving light quarks have to 
be taken into account using the relevant color- 
octet OMEs. Furthermore, contributions from 
Xcj and ip' decays have to be included, whereas 



those from B meson decays are suppressed by the 
small branching fraction. 

The differential cross section of e + e~ — » 
e + e~ H + X, where H denotes a generic charmo- 
nium state, can be written as 



~H + X) 



x 



a.b.d 



x / dx b f bh (x b ,fj, f )y2(0 H [n]) 



x da(ab — ► cc[n] + d), (1) 

where f 1 / e { x ±) is the equivalent number of 
transverse photons radiated by the initial-state 
positrons and electrons [ |l3|, f a /j(x a , [if) are the 
PDFs of the photon, (O^n]) are the OMEs of 
the H meson, da(ab — > cc[n] + d) are the dif- 
ferential partonic cross sections, the integrals are 
over the longitudinal-momentum fractions of the 
emitted particles w.r.t. the emitting ones, and it is 
summed over a, 6 = j,g,q,q and d = g,q,q, with 
q = u,d, s. To leading order in v, we need to in- 
clude the cc Fock states n = 3 s[ 1] , ^ , 3 sf ] , 3 P l f ] 

if H = J/ip,tp' and n = 3 P [ j\ 3 sf ] if H = 
Xcj, where J = 0,1,2. With the definition 
/ 7 / 7 (x 7 , p,f) = 5(1 — Xy), Eq. ([!]) accommodates 
the direct, single-resolved, and double-resolved 
channels. 

In our numerical analysis, we use m c — (1.5 ± 
0.1) GeV, a = 1/137.036, and the LO formula 
for a{ nf \p) with rif — 3 active quark flavors. As 
for the photon PDFs, we use the LO set from 
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Figure 2. Transverse- momentum distribution da/dpr of 77 — * J /if) + X via bremsstrahlung at CERN 
LEP2. The sum of the LO contributions for X = j and X = 7 is compared with the 2^3 part of the 
NLO contribution for dijet invariant mass Sjj > M 2 and M 2 /20. For comparison, also the 3 s[ 8 ' -channel 
contributions to the latter arc shown. 



Ghick, Reya, and Schienbein (GRS) [ |l4j], which 
is the only available one that is implemented in 
the fixed-flavor-number scheme, with nf = 3. 
We choose the renormalization and factorization 
scales to be \i = ^tot and (if = £firiT, respec- 
tively, where m<r = */ M 2 + pip is the transverse 
mass of the J/ij) meson, and independently vary 
the scale parameters £ M and £/ between 1/2 and 
2 about the default value 1. As for the J/tp, XcJ, 
and ip' OMEs, we adopt the set determined in 
Ref. [ [[5) by fitting the Tevatron data using the 
LO proton PDFs from Martin, Roberts, Stirling, 
and Thorne (MRST98LO) [|l6j as our default and 
the one referring to the LO proton PDFs from the 



CTEQ Collaboration (CTEQ5L) [ for compar- 
ison (see Table I in Ref. In the first (second) 
case, we employ Aqq D = 204 MeV (224 MeV), 

which corresponds to A^ D = 174 MeV [ |l6) 
(192 MeV [ 0), so as to conform with the fit [ 
[L5| . Incidentally, the GRS photon PDFs are also 
implemented with A^ D = 204 MeV [[H). In the 
cases V = J/i>,ip', the fit results for (0* fS^ 1 ]) 
and (O^PPq 8 ']) are strongly correlated, and one 
is only sensitive to the linear combination 



Mf = (O^sf]) + — (0*[ 3 P' S| ]), 
mi 



■yip [3p[8]i 



(2) 
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with an appropriate value of r. Since Eq. ([!]) 
is sensitive to a different linear combination 
of (O^Sf 1 ]) and (C^pPff 1 ]} than appears in 

Eq. (|), we write (O^Sf 1 ]) = nMf and 
— k) (m^/r) M,r and vary k be- 
tween and 1 about the default value 1/2. The J- 
dependent OMEs (C^ [ 3 Pj 8] ]}, (0*" [ a Pj 1] ]), and 

(P XcJ psj 8 ']) satisfy multiplicity relations, which 
follow to leading order in v from heavy-quark 
spin symmetry. In order to estimate the theo- 
retical uncertainties in our predictions, we vary 
the unphysical parameters £ M , £/, and k as indi- 
cated above, take into account the experimental 
errors on m c , the decay branching fractions, and 
the default OMEs, and switch from our default 
OME set to the CTEQ5L one, properly adjust- 



.(3) 



We then combine the individual shifts 



m S A QCD- 

in quadrature, allowing for the upper and lower 
half-errors to be different. 

In Fig. H w e confront the p\ distribution of 
e+e~ -> e + e~J/ip + X measured by DELPHI [ 
||] with our NRQCD and CSM predictions. The 
solid lines and shaded bands represent the central 
results, evaluated with our default settings, and 
their uncertainties, respectively. We observe that 
the DELPHI data clearly favors the NRQCD pre- 
diction, while it significantly overshoots the CSM 
one [@. 

4. CONCLUSION 

In summary, the production of J/tp mesons in 
77 collisions allows for a test of factorization in 
color-octet processes as predicted by NRQCD and 
observed at the Tevatron. We have calculated 
the contributions from direct and resolved pho- 
tons to J/ip mesons produced directly or from 
Xcj and tj;' decays. Our NRQCD predictions are 
nicely confirmed by recent data from the DEL- 
PHI collaboration at CERN LEP2, whereas the 
CSM prediction is clearly disfavored. 

ACKNOWLEDGMENT 

The author thanks the organizers of the RAD- 
COR/LL 2002 conference for the kind invitation 
and B.A. Kniehl, L. Mihaila, and M. Steinhauser 
for a fruitful collaboration. 



REFERENCES 

[I] H. Fritzsch, Phys. Lett. B 67 (1977) 217. 
[2] F. Halzen, Phys. Lett. B 69 (1977) 105. 

[3] M. Gliick and E. Reya, Phys. Lett. B 79 
(1978) 453. 

[4] P. Hoycr and S. Peigne, Phys. Rev. D 59 

(1999) 034011. 

[5] N. Marchal, S. Peigne and P. Hoyer, Phys. 

Rev. D 62 (2000) 114001. 
[6] E. L. Berger and D. L. Jones, Phys. Rev. D 

23 (1981) 1521. 
[7] R. Baier and R. Riickl, Phys. Lett. B 102 

(1981) 364. 

[8] W. E. Caswell and G. P. Lepage, Phys. Lett. 

B 167 (1986) 437. 
[9] G. T. Bodwin, E. Braaten and G. P. Lepage, 

Phys. Rev. D 51 (1995) 1125 [Erratum-ibid. 

D 55 (1997) 5853]. 
[10] M. Klasen, B. A. Kniehl, L. Mihaila and 

M. Steinhauser, Nucl. Phys. B 609 (2001) 

518. 

[II] M. Kramer, J. Zunft, J. Steegborn and 
P. M. Zerwas, Phys. Lett. B 348 (1995) 657. 

[12] S. Todorova-Nova, in Proceedings of the 
XXXI International Symposium on Multipar- 
ticle Dynamics, Datong, China, 1-7 Septem- 
ber 2001 (World Scientific, Singapore, to ap- 
pear) ]hep-ph/0112050C ; M. Chapkin, talk 
presented at 7th International Workshop on 
Meson Production, Properties and Interac- 
tion (Meson 2002), Krakow, Poland, 24-28 
May 2002; private communication. 

[13] S. Frixione, M. L. Mangano, P. Nason and 
G. Ridolfi, Phys. Lett. B 319 (1993) 339. 

[14] M. Gliick, E. Reya and I. Schienbein, Phys. 
Rev. D 60 (1999) 054019 [Erratum-ibid. D 62 

(2000) 019902]. 

[15] E. Braaten, B. A. Kniehl and J. Lee, Phys. 

Rev. D 62 (2000) 094005. 
[16] A. D. Martin, R. G. Roberts, W. J. Stirling 

and R. S. Thome, Eur. Phys. J. C 4 (1998) 

463. 

[17] H. L. Lai et al. [CTEQ Collaboration], Eur. 

Phys. J. C 12 (2000) 375. 
[18] M. Klasen, B. A. Kniehl, L. N. Mihaila and 

M. Steinhauser, Phys. Rev. Lett. 89 (2002) 

032001. 



5 



NRQCD 

3 p [8] 

3p J [l] 

lo J [8] 
^0 



e + e -> e + e J/\\f X at LEP2 



DELPHI prelim. 
MRST98 fit 

CTEQ5 fit 
VS = 197 GeV 

- 2 <y J/¥ < 2 




Pf (GeVO 



Figure 3. The cross section da jdp\ of e + e~ — > e + e~ J/i/j + X measured by DELPHI [ 12 as a function 
of is compared with the theoretical predictions of NRQCD and the CSM. The solid and dashed lines 
represent the central predictions obtained with the ME sets referring to the MRST98LO [ |l6| (default) 
and CTEQ5L [ |l^| PDFs, respectively, while the shaded bands indicate the theoretical uncertainties on 
the default predictions. The arrows indicate the NRQCD prediction for px = and its 3 P^\ ^(j 8 ', 3 S^\ 
and 

3p[8] 

components. 



